Controlling nuclear spin exchange via optical Feshbach resonances in ^'^^Yb 
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Nuclear spin exchange occurs in ultracold collisions of fermionic alkaline-earth-like atoms due to 
a difference between s- and p-wave phase shifts. We study the use of an optical Feshbach resonance, 
excited on the ^5*0 ^Pi intercombination line of ^^^Yb, to affect a large modification of the s- 
wave scattering phase shift, and thereby optically mediate nuclear exchange forces. We perform 
a full multichannel calculation of the photoassociation resonances and wave functions and from 
these calculate the real and imaginary parts of the scattering length. As a figure of merit of this 
interaction, we estimate the fidelity to implement a VSWAP entangling quantum logic gate for two 
atoms trapped in the same well of an optical lattice. For moderate parameters one can achieve a 
gate fidelity of ^ 95% in a time of ^ 50/is. 
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Interest in the use of ultracold alkaline-earth-like atoms 
for applications in quantum information processing has 
grown in recent years due to some unique features of 
their electronic structure and experimental advances in 
quantum control for optical-clock technology [ ] . Because 
the ground state is a closed shell state, for fermionic 
species, quantum coherence can be stored in the nuclear 
spins that are isolated from hyperfine coupling and very 
weakly coupled to perturbing noisy magnetic fields. The 
nuclear spins, nonetheless, can mediate strong electronic 
interactions between atoms due the quantum statistics 
of identical particles [ ]. In addition, the existence of 
meta-stable excited states provides further avenues 
for control by transferring quantum coherence between 
nuclear and electronic degrees of freedom. These various 
tools have been assembled in a number of proposals for 
quantum control, including cooling of atomic vibration 
without decohering nuclear spin coherence [ ], architec- 
tures for quantum computers [ ], and quantum simula- 
tions of exotic forms of quantum magnetism [ ]. 

An important additional tool for quantum manipula- 
tions is the capability to control the interaction strength 
between atoms through Feshbach resonances [ ]. This 
has been essential in the study of quantum many-body 
phenomena in dilute ultracold gases [7, 8], enabling tune- 
able interactions for applications including the study of 
quantum phase transitions such as the BEC-BCS cross- 
over [ ], and the production of ultracold molecules [ ]. 
In systems of cold alkali atoms, Feshbach resonances oc- 
cur in the ground electronic manifold between open and 
closed hyperfine channels, and are typically tuned with 
magnetic fields, though optical tuning has also recently 
been demonstrated [ ] . The lack of hyperfine structure 
in the ground state of the alkaline-earths removes this 
possibility, but the existence of the meta-stable ^P ex- 
cited states makes these atoms attractive for implemen- 
tation of optical Feshbach resonances (OFR) [12, 13, 14], 
induced through laser coupling of the scattering atoms 
to an electronic-excited-state bound molecule. In con- 




FIG. 1: (Color online) Schematic of the OFR, leading to nu- 
clear spin exchange. Optical excitation of atoms colliding in 
an s-wave is possible only for the nuclear spin singlet (a) . Nu- 
clear spin triplet p-wave collisions (b) are suppressed at low 
energies and because of selection rules, since they connect to 
different molecular bound states (see text). 



trast to alkali atoms where observation of an OFR is 
accompanied by substantial losses due to rapid sponta- 
neous emission [ , ^"^J, the narrow intercombination 
lines are ideal for use in an OFR. Optical effects in colli- 
sions of alkaline-earth-like atoms have been demonstrated 
in recent experiments, including photoassociation spec- 
troscopy in ^'^^Yb [ ] and ^^Sr [19] and OFR of several 
bosonic isotopes with zero nuclear spin, including ^^^Yb, 
Yb [ ]and^8Sr[ ]. 

In this letter we study OFRs as a tool for optically 
controlling entangling interactions between nuclear spins 
of ^^^Yb atoms, a fermionic group-II-like element with a 
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spin 1 = 1/2 nucleus. This control can be accomplished 
through the laser modification of s-wave scattering in- 
teraction. Because of the exchange symmetry of iden- 
tical fermions, s-wave collisions (dominant at ultracold 
temperatures) are allowed only for the anti-symmetric 
nuclear-spin singlet pairing. A relative phase between 
nuclear-spin singlet and triplet of 7r/2 arising in the col- 
lision generates an entangling ^/SWAP^ studied previ- 
ously by Hayes et at. [ ], and observed by Anderlini et 
al. [22] in an equivalent scheme for hyperfine states in 
alkali atoms. In this paper we show how the state of the 
nuclear spins allows or forbids the optical coupling of the 
colliding atoms to an excited molecular bound state, and 
thereby mediates the presence of an OFR, which strongly 
affects the exchange force (see Fig. 1). 

We consider OFRs based on laser excitation near the 
strongest of the intercombination lines, ^5*0 ^ ^Pi. To 
model the Feshbach resonance, we require a good under- 
standing of both the ground-state scattering wave func- 
tion and the bound-excited- molecular wave functions. 
For alkaline-earth-like atoms, there is a single electronic 
^5*0 + ^5*0 ground state channel. For scattering states at 
low energies, the bulk of the probability amplitude of the 
wave functions lies away from the chemical binding re- 
gion, and we need not model the ^5*0 + ^Sq potential 
in detail. Instead, we make use of the analytic solu- 
tions to a van der Waals potential found by Gao 
[23]. A superposition of the regular and irregular solu- 
tions, matched by the boundary conditions to achieve the 
correct asymptotic solution with the appropriate known 
scattering length [ ], yields an excellent approximation 
to the exact wave function for separations larger than 
- lOA. 

The excited states are much more complicated, giv- 
ing rise to a multitude of different coupled potentials, 
and thus requiring a full multichannel description. The 
Hamiltonian for the diatomic molecule that asymptotes 
to unbound ^6*0 + ^Pi atoms has the form 

pI 

^ = ^ + Hrot + HbO + HhF- (1) 

The first two terms represent the radial and angular ki- 
netic energy of the nuclei. Hbo represents the elec- 
tronic Born-Oppenheimer potentials that asymptote to 
the ^5*0 + ^Pi scattering channel, taken in "Hund's-case 
(c)" where the spin-orbit interaction is the dominant 
effect. This is modeled through a Lennard- Jones plus 
dipole-dipole potential, 

^^7.(r) = ^-^-a^, (2) 

with parameters determined by fits to experiments as 
Ce = 2810 au, C12 = 1.862 x 10^ au and Cl = = 
0.09695 au. Hhf is the hyperfine interaction in the ^Pi 
orbital. We neglect the small magnetic dipole-dipole in- 



teractions between nuclei. We do not include an external 
magnetic field in this calculation. 

To determine the eigenstates, we perform a multichan- 
nel diagonalization. Note that for the given Hamiltonian 
there are only a few exact good quantum numbers: the 
total angular momentum quantum number T, its projec- 
tion on the space- fixed quantization axis Mt, and par- 
ity due to inversion of all particles p. Other quantum 
numbers are approximate, depending on the dominant 
forces. Three different channel bases defined by differ- 
ent couplings of the angular momenta of the atoms are 
useful in analyzing the eigenstates [25]. The Hbo is di- 
agonal in a Hund's case (c) basis, extended to include 
nuclear spin, I7) = \JQIl^{T^ Mt^p))^ where J^I^F are 
the magnitudes of the electron, nuclear, and total (ex- 
cluding rotation) angular momentum, with projections 
on the internuclear axis 1], ^, <l>, respectively. In contrast, 
the rotational and hyperfine Hamiltonians are diagonal 
in an extended Hunds case-(e) basis which is separa- 
ble between nuclear rotation and other degrees of free- 
dom, |e) = \hf2FR{T,MT,p)). Here h = 1/2 and 
/2 = 1/2,3/2 are the possible individual atomic an- 
gular momenta for the ground and excited state, re- 
spectively. A third basis of importance is the prod- 
uct basis of internal atomic states and nuclear rotations, 
k) = |/i,mi,/2,m2,P(T,MT,p)). 

Selection rules divide the excited-state channels into 
two classes: those that are optically accessible by ground- 
states colliding in an s-wave from those accessible in a 
p-wave collision. Because of the correlation between nu- 
clear spin pairing and collisional partial wave, these se- 
lection rules play an essential role in optical control of 
nuclear exchange forces. The dipole-allowed transition 
requires that the ground and excited states have oppo- 
site parity and that total angular momentum changes by 
AT = 0,1 (T = 0^0 is forbidden). By exchange sym- 
metry, in the ground state, the two- atom nuclear spin sin- 
glet (/ = 0) is associated with s-wave collisions {R = 0) 
and therefore has even parity p = and total angular mo- 
mentum T = 0. Thus the excited state must have p = —1 
and T = 1. Taking the selection rules into account, there 
are 5 excited channels, each 3-fold degenerate in the ab- 
sence of an external magnetic field, which can be opti- 
cally excited for an s-wave collision in the ground state. 
In contrast, the two- atom nuclear spin triplet (/ = 1) 
is associated with p-wave collisions {R = 1) and there- 
fore we have odd parity and possible total angular mo- 
menta T = 0, 1, 2 in the ground state. By the selections 
rules, the electric-dipole connected excited states must 
have T = 0, 1, 2, 3, and even parity. In total, the p-wave 
collisions are optically coupled 19, (2T+l)-fold degener- 
ate channels in the absence of a magnetic field (a total 
of 89 channels). 

A full multichannel diagonalization of the Hamiltonian 
based on the DVR method [ ] yields the bound states 
representing photoassociation resonances that dissociate 
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to ^S'o + = 3/2). The s-wave-accessible resonances 

are shown Table 1. These are in good agreement with 
the experimental observations of Enomoto et al. [ ], 
though some states were not reported in that experi- 
ment, either because the binding energy is comparable 
to the linewidth or because these state are too tightly 
bound and thus have small Franck- Condon factors. We 
also note the existence of potential curves that support 
so-called purely long-range (PLR) states, observed by 
Enomoto et al. for p- waves. These potentials arise from 
avoided crossings due to hyperfine mixing of different 
Hund's case-(c) potentials. Because of the weakness of 
this interaction, the avoided crossings occur at very large 
internuclear separations and the resulting potential is ex- 
tremely shallow (a 59.2 MHz deep PLR potential acces- 
sible through laser excitation of an s-wave collision and a 
717 MHz deep PLR potential for p-waves). Nonetheless, 
they can support bound molecular states. 

With the so determined wave functions and eigenener- 
gies, we calculate the properties of our OFR. The optical 
modification of the s-wave scattering length can be ex- 
pressed as a(/. A) = ai,g + aopt{I,^) - ibopt{I,^) [12], 
where a^^ is the background scattering length, 

is the optical contribution to the elastic scattering length, 
depending on the intensity / of the radiation and the 
detuning A, and 

A) = ^ (^^, ^ Ir^uUijf/^) ' 

is the imaginary part of the scattering length, accounting 
for two-body scattering loss. The molecular spontaneous 
decay rate, Fm, is calculated from the multichannel wave 
function. The linewidth is broadened by stimulated emis- 
sion at the rate [ ], 

rstUi) = I hTlUotfFc, (5) 

Here F/^/27r = 182 kHz is the atomic natural linewidth 
for the ^Sq ^ ^Pi intercombination transition in ^^^Yb, 
Isat = {27r^hrAc)/{3X^) = 0.13 mW/cm^ is the atomic 
saturation intensity. The Franck- Condon factor that 
measures the overlap between the spatial wave function of 
the (energy normalized) scattering ground electron state 
at energy E and the (unit normalized) bound excited 
state with the appropriate rotational modification [ ] is 
given by 

fFcfrot = \{i^e\d-eL\ijg)f /2dl, (6) 

where dA is the atomic dipole moment for this transition 
and e^;, is the laser polarization. Equation (6) is calcu- 
lated using the multichannel wavefunctions expanded in 
the basis |7r), as discussed above and in [^u]. 



The key parameter that determines the strength of the 
resonance is the "optical length" defined by lopt{I) = 
^ stim{I) / {'^K^ m) ^ where kr is the wave number for the 
relative motion of the reduced mass /i. Table I gives the 
optical length for the s-wave-accessible photoassociation 
resonances near dissociation, calculated at a laser inten- 
sity of 1 W/cm^. Given these values, we can find the 
optically modified scattering length as a function of laser 
intensity and detuning. Since the bare s-wave scattering 
length of ^^-"^Yb, a^^ = — 0.15ao, is essentially zero, an 
optical Feshbach resonance will increase |a| by orders of 
magnitude. 



The utility of the OFR for coherent control of nuclear 
spin exchange depends on low loss and decoherence. As 
a figure of merit, we consider the example of implement- 
ing a V SWAP two-qubit entangling unitary between the 
spin-1/2 nuclei trapped in one site of an optical lattice, 
as in [22]. Such a gate occurs if the relative phase shift 
for singlet vs. triplet nuclear spin states is 7r/2. In the 
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24.7 


245 


26887 


34.9* 




89 


7382 


49.9 
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89.1 


86.9 


248 


5953 


152.3 


149.9 


251 


3198 


250.2 


247.0 


256 


1495 


396.5 


391.7 


261 


840 


462.1 




161 


13 


597.3 




101 


235 


612.5 


604.3 


241 


390 



TABLE I: S-wave accessible excited molecular bound states 
and resonance properties. ^^^^°^y ^^^e the binding energies 
found from the multi-channel calculation solutions of the 
Schrodinger equation with Hamiltonian Eq. (1). PLR states 
are marked with *. Binding energies are denoted in frequency 
units relative to the atomic ^ Sq ^Pi{f = 3/2) transition. 
For comparison, E^^-^ are the photoassociation resonances ob- 
served by Enomoto et al. with an experimental uncertainty 
of ±2MHz. The agreement between experiment and theory is 
good, with a small systematic shift that may be attributed to 
inaccuracies in the Hund's case-(c) potentials, and/or to sys- 
tematic shifts in the experimental conditions (e.g. light-shift 
from the confining dipole potential). Fm is the molecular nat- 
ural linewidth and lopt is the optical length at an intensity of 
1 W/cm^ 
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^^^^ (MHz) 


212.0* 


212.4 


233.8* 


234.0 


258.2 


256.9 


270.4 


268.3 


278.9 


276.8 


355.3* 


355.4 


383.2* 


383.4 


415.5 


416.1 





E^^ (MHz) 


431.9 


432.0 


443.2 


442.5 


646.3 


646.2 


667.7 


667.2 


682.4 


681.8 


976.0 


976.2 


1003.5 


1002.1 


1022.1 


1021.4 



TABLE H: P-wave accessible excited molecular bound states. 
EI^^° are the binding energies found from the multi-channel 
calculation as in Table I. Only the lines which were also ex- 
perimentally observed by Enomoto et al. at a temperature 
of 25/iK are shown. The experimental uncertainty is zb2MHz 
and ±1 MHz for the PLR states denoted by *. 



perturbative regime and neglecting the small background 
scattering length, the collisional s-wave phase shift after 
interaction for time T is = {Aiihnaopt/ l~^)T ^ where n is 
the atomic density. The loss/decoherence rate is given by 
2Kn = Snhnbopt/ fJ^' The fidelity of the gate based only 
on this loss is then 



-2KnT 



-(f)(^) 



(7) 



At the collision energies we consider, p-wave phase shifts 
are negligible, even in the presence of laser excitation. 

To get an estimate of how well the OFR induces nuclear 
spin exchange, we must balance a variety of constraints. 
High fidelity at moderate intensities requires large de- 
tuning from molecular resonance, e.g., a fidelity of 95% 
is achieved when A ^ SOTm- Our model, however, as- 
sumes sufficiently small detuning so that a single molec- 
ular excited state contributes to the resonance. In addi- 
tion, to ensure a reasonably fast interaction, the applied 
intensity must be sufficiently large so that the optical 
scattering length is large, which will power-broaden the 
atomic resonance Va \/l + I /Isat^A- For all of these 
reasons, we consider as an example the photoassociation 
resonance bound by 396.5 MHz, with a natural linewidth 
of Tm = 261 kHz. This bound state is still relatively close 
to dissociation, and thus the molecule is long-range, with 
an outer turning point at 130 ao- Nonetheless, the density 
of states is sufficiently sparse that one can detune many 
linewidths from that resonance while still neglecting cou- 
pling to the next higher molecular state, which is bound 
by 250.2 MHz and has about twice the Franck-Condon 
factor. A detuning from the molecular resonance of 
A = -SOFm ~ -7.8MHz fits this constraint. The inten- 
sity is chosen to broaden the resonance so as to increase 
the scattering length by fixing TsUm = '^krlopt^M — —A, 
or /opt ^ 10^ cLq. From Table I, this is achieved at an in- 
tensity of / = 12.7 W/cm^, whereby the atomic linewidth 
is power broadened to V a 57 MHz, which is still nar- 
row compared to a detuning of 400 MHz from dissoci- 
ation. When — A = V sUm ^ Fm, the optical scatter- 
ing length is aopt = —281 ao and the loss coefficient is 



K ^ 1.15 X 10~^^ cm^/s. With this large magnitude 
scattering length and low loss, given two ^^^Yb atoms 
in a lattice site analogous to the experiments at NIST 
[22] with a density n = 2.4 x lO^cm"^, the time of the 
VSWAP gate isT ^ 50/is and the gate fidelity is ^ 95%. 

In principle, higher fidelity, larger scattering lengths, 
and shorter gate times are possible by using a higher in- 
tensity and larger detuning, though a proper treatment 
will require the modelling of excitation to multiple molec- 
ular bound states and line broadening. Even with these 
modest parameters, we see that OFRs have great poten- 
tial for control of nuclear-spin exchange and strong en- 
tangling interactions. The combination of this tool with 
recent advances in loading optical lattices via superfluid 
to mott insulator phase transition [27], the ability to op- 
tically manipulate nuclear spin coherence [ ] and re-cool 
atoms without decohering nuclear spins [ ], and propos- 
als for quantum logic [ , ] make this system attractive 
for new applications in quantum information processing. 
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